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ABSTRACT 



Context. Observation of star occultations is a powerful tool to determine shapes and sizes of asteroids. This is key information necessary for 
studying the evolution of the asteroid belt and to calibrate indirect methods of size determination, such as the models used to analyze thermal 
infrared observations. Up to now, the observation of asteroid occultations is an activity essentially secured by amateur astronomers equipped 
with small, portable equipments. However, the accuracy of the available ephemeris prevents accurate predictions of the occultation events for 
objects smaller than ~100 km. 

Aims. We investigate current limits in predictability and observability of asteroid occultations, and we study their possible evolution in the 
future, when high accuracy asteroid orbits and star positions (such as those expected from the mission Gaia of the European Space Agency) 
will be available. 

Methods. We use a simple model for asteroid ephemeris uncertainties and numerical algorithms for estimating the limits imposed by the 
instruments, assuming realistic CCD performances and asteroid size distribution, to estimate the expected occultation rate under different 
conditions. 

Results. We show that high accuracy ephemerides which will be available in the future will extend toward much smaller asteroids the possibil- 
ity of observing asteroid occultations, greatly increasing the number of events and objects involved. A complete set of size measurements down 
to ~10 km main belt asteroids could be obtained in a few years, provided that a small network of ground-based lm telescopes are devoted to 
occultation studies. 
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ctf 1. Introduction 

An accurate knowledge of asteroid sizes and size distribution 
is a fundamental prerequisite for studying the origin and the 
collisional evolution of these bodies, including the formation 
and the dispersion of dynamical families. It also has important 
implications for our understating of the mechanisms that led 
to the formation of the terrestrial planets of our Solar System. 
Furthermore, knowledge of the size of an asteroid allows its 
bulk density to be inferred when the mass of the body can be in- 
dependently determined, for instance from the orbital period of 
an asteroid satellite. Information about the density of asteroids 
is crucial to study their interiors, and to constrain the collisional 
processes leading to the formation and evolution of the differ- 
ent populations of these minor bodies (e.g. IBritt et al. 20 02 ). 
Moreover, if the size of an asteroid is known, the albedo of 
the body can be derived using Eq. Q](e.g. Bow ell et al. 19 89 ), 



where pv is the geometric visible albedo and H is the absolute 
magnitude of the asteroid. 



D(fen) = ^ x 10- H/5 



(1) 
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Note that H is a quantity that can be easily obtained from mea- 
surements of the apparent brightness of the object, for exam- 
ple by means of photometric observations in the visible light. 
The distribution of albedos and its correlation with the taxo- 
nomic types and the orbital elements of asteroids is a key piece 
of data in understanding their nature and their origin, and can 
help to constrain the mineralogical composition of their sur- 
faces. Albedos are also crucial for developing and improving 
reliable and truly debiased asteroid population models (e.g. 
ITedesco et al. 2005IIStuart and Binzel 20041 ). 

Asteroid flybys of space probes have provided not only the 
best "direct" size determination, but also a wealth of informa- 
tion on shapes and surface properties for a very limited num- 
ber of these bodies. However, a statistically significant sample 
of asteroid sizes can be obtained only by remote observations. 
Up to now, mainly indirect methods have been used to deter- 
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mine asteroid diameters: the vast majority of asteroid sizes and 
albedos have been derived from observations obtained in the 
thermal infrared (~ 12 - 100 pmi) by the Infrared Astronomical 
Satellite (IRAS; |Tedesco 19921 ITedesco et al. 20021 . However, 
it is known that IRAS diameters of asteroids are not exempt of 
systematic errors. They were derived by means of the Standard 
Thermal Model (STM of |Lebofsky and Spencer"T989| , whose 
results are dependent on the value of a model parameter called 
the beaming factor and indicated with 77. The value of 77 was 
calibrated by forcing the STM to give the correct diameters 
(as derived from occultation observations) of the asteroids 1 
Ceres and 2 Pallas (Lebofsk y et al. 1986| l. Several studies have 
highlighted, though, that a value of r\ appropriate for large 
main belt asteroids, might lead to an underestimation of the 
true size when applied to small asteroids whose surface ther- 
mal characteristics are likely to be different from those of 
larger bodies (e.g. |Spencer et al. 1989] IDelbo & Harris 20021 
Harris and Lag erros 2002][ 

Calibration of asteroid radiometric diameters was possible 
because, in the past half century, the sizes of the largest aster- 
oids were also derived by measuring the duration of star occul- 
tations. This approach has traditionally been flawed by the poor 
accuracy of the prediction of the asteroid shadow path on Earth, 
thus keeping the success rate of asteroid occultation very low. 
The availability of precise star astrometry obtained by the mis- 
sion Hipparcos of the European Space Agency (ESA), has im- 
proved the situation in the last ~ 10 years (Dunh am et al. 20 02 ), 
strongly increasing the number of successfully measured aster- 
oid occultation chords. Electronic equipments (fast video cam- 
eras) have also increased the accuracy of occultation timings, 
previously obtained mainly by visual observation. However, as- 
teroid occultation observations have remained mainly an activ- 
ity reserved to observers equipped with small, movable tele- 
scopes and equipment, ready to accept a rather high rate of 
"miss". 

In this study, after having reviewed the present state of as- 
teroid occultation studies, we will focus on a very likely evolu- 
tion of this observational approach, thanks to the future avail- 
ability of asteroid orbits and star positions whose uncertainties 
will be reduced by orders of magnitude. In this respect, the 
astrometric ESA mission Gaia (Perryman et al. 2001 ) will rep- 



resent a major milestone that will completely change the ap- 
proach to asteroid occultation studies, opening the way to the 
possibility of obtaining a more complete mapping of asteroid 
sizes and albedos. In this paper we address this forthcoming 
improvement, exploring the path that the approach to asteroid 
occultations will take in the future and investigating the limits 
that future observation will be able to approach. 

Note that indirect methods of asteroid size and albedo de- 
termination, such as those based on polarimetry and thermal in- 
frared observations, will greatly benefit from occultation stud- 
ies performed on a larger sample of asteroids: it will be possi- 
ble, for example, to calibrate models used to analyze thermal 
infrared and visual polarimetric measurements on a larger and 
more accurate sample of calibration asteroid diameters, with 
object sizes extending down to some tens of kilometers. 

Since better performing telescopes will be needed to push 
size determination toward smaller and thus fainter asteroids, in 



this paper, we also study the observing rate of asteroid occulta- 
tions expected for a fixed observing station using a one-meter- 
class telescope equipped with a cooled CCD camera capable of 
a high frame rate. 

After having assessed the present state of occultation stud- 
ies (Sect. 1) we analyze the theoretical occultation rate that is 
expected for different asteroid sizes (Sect. 2). We then study 
the role of instrument detection capabilities (Sect. 3) and that 
of ephemeris accuracy (Section 4). The implications of high- 
accuracy ephemeris are then discussed (Section 5). 

2. The observation of asteroid occultation: where 
we are now 

2.1. Observational campaigns 

The first occultation ever recorded is that of 3 Juno, in 1958. 
The third one (433 Eros), observed in 1975, was also the first 
to be observed by multiple sites, allowing the determination of 
different "chords" of the asteroid (Millis and Dunham 1989). 
Only 17 positive occultation observations were obtained by the 
end of 1980. 

Coordinated campaigns and last-minute small-field astrom- 
etry, aimed both at improving the occultation path and at dif- 
fusing the relevant predictions among observer networks, have 
provided a moderate increase in the rate of successfully mea- 
sured asteroid occultations events. Up to the early 1990s the 
probability of observing an occultation chord for a given event 
was discouraging. However, consistent efforts were rewarded 
with a total of 60 positive events by 1985, growing to 101 in 
1990, and increasing more rapidly later, thanks to the avail- 
ability of more precise star astrometry. After the publication 



of the Hipparcos and Tycho-2 star catalogues (H0g et al. 2000) 
the number of successfully recorded occultation events doubled 
in the decade 1991-2000 (see below for an evaluation of the 
current prediction accuracy). At the end of 2004 they summed 
up to 658, and today they continue to grow at a rate of about 
-100 per year dDunham & Herald 2005b . 

Despite the apparently large number of observed events, the 
extremely limited number of chords for most of the observed 
events results in diameters that remain poorly constrained. 
Also, the poor knowledge of shapes, rotational poles, and ro- 
tational phases at the epoch of the events often makes the in- 
terpretation of observations difficult. In practice, the use of oc- 
cultation diameters found in current dataset^] is not straight- 
forward. As an example, IShevchenko and Tedesco 20 06 were 
able to use just 57 occultation diameters for albedo determi- 
nations, among which only 18 have an estimated accuracy 
< 5%.fShevchen ko and Ted esco 2006 compared their "occulta- 
tion albedos" against the albedos derived from IRAS observa- 
tions in the thermal infrared, and against the presently available 
polarimetric albedos. However, because their data set contains 
occultation sizes only for asteroids larger than ~100 km, they 
were not able to explore any size-dependence of systematic er- 
rors of polarimetric and thermal infrared albedos. 



1 Interested readers should refere to the PDS asteroid occultation 
dataset (AOD), available from the NASA Planetary Data Systems 
Small Bodies Node: http://www.psi.edu/pds/archive/occ.html 
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Asteroid occultation observations - strongly encouraged by 
professional astronomers - are supported and secured by or- 
ganizations grouping mainly volunteering amateurs, such as 
the International Occultation Timing Association (IOTA) or 
the European Asteroid Occultation Network (EAON), strictly 
cooperating with - and including - professional astronomers. 
Today, for stars currently included in the predictions of as- 
teroid occultations (usually brighter than V=12 and listed in 
UCAC or Hipparcos catalogues) the main source of prediction 
uncertainties resides in the accuracy of asteroid ephemerides 
(as discussed in the next section). Current uncertainties in the 
ephemerides of asteroids are rarely below ~0.5 arcsec (ex- 
cept for the largest bodies), corresponding to a displacement 
of ~350 km on the Bessel plane for an object at 1 AU from 
the Earth. The corresponding uncertainty on Earth has similar 
values, thus requiring a large number of appropriately spaced 
observers to increase the success probability of actually observ- 
ing the occultation. Furthermore, in the coming years a slight 
degradation in the accuracy of occultation prediction can be 
expected, as a consequence of the uncertainty on the proper 
motions of stars. 

The Euraster network data^l show that in the year 2005, 
~15% of the occultations observed under a clear sky was suc- 
cessful. This simple evaluation, however, does not include any 
information on the real position of the observers with respect 
to the predicted path on Earth, nor a distinction based on the 
presence of other degrading factors, such as a small flux drop, 
ephemeris uncertainty, short duration, and bad atmospheric 
conditions. It would thus be difficult and time-consuming to es- 
tablish a more precise "average" success rate, since it depends 
on details that would need a case-by-case evaluation. 

2.2. Current accuracy of asteroid ephemeris 

The uncertainty in asteroid ephemeris can be evaluated by 
using information provided by the Bowell's Asteroid Orbit 
Database (ftp://ftp.lowell.edu/pub/elgb/astorb.html). The aster- 
oid ephemeris uncertainty in this source is calculated by means 
of the procedure described by Muinonen and Bowell (1993), 
and it is expressed in terms of a small set of parameters, namely 
the CEU (Current Ephemeris Uncertainty), the PEU (Peak 
Ephemeris Uncertainty, and its date of occurrence), and the 
PEU over the following 10 years. All quantities are referred to 
the constantly updated epoch of the osculating orbital elements. 
The method adopted for the calculation of the CEU is consid- 
ered to give meaningful estimates of the current ephemerides 
uncertainty for multi-opposition Main Belt Asteroids (MBAs) 
but, on the other hand, it is assumed to be unreliable for Earth 
crossers or asteroids with poorly determined orbits. Since we 
focus our study on the observation of occultations events of 
MBAsH we can assume that the CEU is representative of the 



2 http://www.euraster.net/results/index.html 

3 According to Tedesco et al. (2005) we define MBA as those aster- 
oids having orbits between the 4: 1 and 2: 1 mean-motion resonances 
with Jupiter, i.e., those with orbital semimajor axes between 2.064 
and 3.278 AU and with modest inclinations (<25°) and eccentricities 
(<0.3). 



uncertainty on the position at the epoch for which the osculat- 
ing orbital elements are calculated. 

Note that the CEU gives the length of the major axis of the 
uncertainty ellipse (usually very slender) representing the set 
of positions on the sky where the object can be found, at a 1- 
sigma confidence level. This axis is usually lying on the line 
of variations, or at a small angle with it. The main effect of 
this orientation will thus be a major uncertainty on the event 
epoch, not on the shadow path in the direction perpendicular to 
the asteroid shadow motion. This latter uncertainty is the dom- 
inating one as far as the observability of the occultation event 
is concerned; in fact, a shift in the predicted epoch can easily 
be compensated by an increased duration of the observation. 
It is thus likely that, by retaining the CEU as an estimate for 
the position uncertainty we are, in general, rather conservative. 
The exception is probably represented by a small number of 
large asteroids, having a smaller and more circular uncertainty 
ellipse, for which the CEU correspond more strictly to the pre- 
diction accuracy (Bowell, private communication). 

Under these assumptions, we can then state that with a CEU 
larger than the apparent angular size of the asteroid, occultation 
events are difficult to predict. We therefore used the ratio of the 
CEU to the apparent angular radius of the asteroid (given 
by the ephemeris at the CEU epoch) as a statistical indicator 
for the current predictability of occultation events. For other 
epochs, of course individual CEU values will be different, but 
the overall statistics remain valid. 

For our computations, the apparent sizes of asteroids 9 
were calculated from their IRAS diameters included in the 
Supplemental IRAS Minor Planet Survey (SIMPS) of Tedesco 
et al. (2002). Note that the SIMPS contains diameters for only 
the first ~2200 numbered asteroids. For the large majority of 
bodies in our sample, diameters are thus not known. For them, 
we estimated sizes from asteroid H magnitudes (which are in- 
cluded in the astorb.dat file) and an assumed geometric visible 
albedo via Eq.Q] However, the use of a constant albedo through 
the entire Main Belt would be an oversimplification: the albedo 
range for MBAs is very large (0.03 - 0.5), and it is known, for 
instance, that an inverse correlation between albedo and helio- 
centric distance exists. We used, therefore, the procedure de- 
scribed by Tedesco et al. (2005) to assign realistic albedos to 
the numbered asteroids not included in the SIMPS database. 

According to this model (and using the numbering as in 
Tedesco et al. 2005), the main belt is divided into three zones: 
i.e. zone 2 for 2.064 < a < 2.50; zone 3 for 2.50 < a < 2.82; 
zone 4+5 for 2.82 < a < 3.278, where a is the orbital semi- 
major axis of the asteroid. To assign an albedo, we first di- 
vide MBAs according to the zone to which they belong (2, 3 or 
4+5) and to their mean apparent visible magnitude at opposi- 
tion, given by V(a, 0) = H + 5 log(a(a - 1)). Either the observed 
or the bias-corrected albedo distribution of Ted esco et al. 2 005 
(Table 3) is assigned to an asteroid, depending on whether 
the value of V(a,0), is < 15.75 or > 15.75, respectively. 
The same table gives the zonal distribution for four different 
classes of albedos: i.e. low (0.020 < pv < 0.089), intermedi- 
ate (0.089 < p v < 0.112), moderate (0.112 < p v < 0.355), 
and high (0.355 < pv < 0.526). We then weight the number of 
asteroids according to the corresponding zonal albedo distribu- 
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Fig. 1. Plot of the ratio CEU/0 versus the object diameter (in km, lower X-axis) for the numbered main belt asteroids known 
and observable at the date of writing: Feb 20, 2007. The solid curve, superimposed on the plot, shows the cumulative number of 
asteroids (upper X-axis) which have a CEU/0 ratio smaller or equal than the corresponding value of the Y-axis on the left-hand 
side of the figure. The curve clearly shows that there are about 80 - 90 asteroids with a CEU/0 value smaller than 1 . We expect 
that Gaia will improve the accuracy of asteroid ephemeris prediction by a factor ~ 100 (dotted horizontal line), increasing by 
almost three orders of magnitudes the number of asteroids for which CEU ~ 0. 



tion and randomly distribute albedos uniformly within each of 
the albedo classes. 

Fig. Q] shows the resulting ratio CEU/0 for all numbered 
MBAs known and observable at the date of writing (astorb.dat 
file of Feb 20, 2007) as a function of their diameter. We as- 
sumed that an asteroid is "observable" (in a broad sense) if it 
has a solar elongation >50". The distribution shows that only 
very few asteroids (less than 90) - mainly bright MBAs with 
diameters larger than 50 km - have a ratio CEU/0 smaller than 
one. This roughly corresponds to the number of asteroids for 
which we can expect today to make reliable predictions of oc- 
cultation events, or stated in another way, their occultations 
have a high probability of being observed inside the predicted 
shadow path. An improvement of occultation predictions is at 
present still possible through last-minute astrometry, aiming to 
constrain the trajectory of the asteroid with respect to the target 
star in the hours immediately preceding the event. 

It is probable that only by the forthcoming all-sky as- 
tronomical surveys, which are capable of providing accurate 



multi-epoch astrometry of asteroids, will we be able to extend 
occultation studies of minor bodies to much smaller sizes, with 
the possibility of deriving occultation diameters for several tens 
of thousand of asteroids. The all-sky survey program that will 
have the best ever astrometric accuracy is the ESA mission 
Gaia, foreseen to be launched at the end of 201 1. 



However, the capability of systematic asteroid occultation 
observing campaigns of providing accurate size determination 
can be judged only if we can estimate the observational effort 
required. Thus, we will calculate in the following section the 
theoretical asteroid occultation rate for a given site on Earth, as 
a function of the asteroid diameter. Later, ephemeris and star 
position uncertainties will be taken into account to estimate the 
success rate of asteroid occultation measurements considering 
instruments with different performances, from those of a 20 cm 
telescope, up to the the case of a one-meter telescope equipped 
with a state-of-the-art CCD camera. 
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3. A numerical model for computing the 
occultation rate 

In order to reproduce a realistic occultation rate for asteroids of 
different sizes, we need to estimate how many stars each aster- 
oid occults in a unit time (e.g. one day). This is equivalent to 
the number of stars contained within the area on the sky swept 
by each asteroid in that time interval, and visible for a generic 
observer on Earth. This area can be easily calculated by tak- 
ing the product between the apparent angular diameter of the 
asteroid, 26, and the distance on the sky between the two po- 
sitions of the body calculated with the difference of one day. 
In order to calculate the density of stars falling within the area 
swept by each asteroid every day, we extracted the stellar den- 
sity from the realistic galactic model used by Crosta & Mignard 
(2006). In this model the number of stars per square degree, as 
a function of their magnitude and their galactic coordinates, is 
expressed in terms of a set of Chebyshev polynomials, whose 
coefficients were determined by a fit to the measured star count. 
Note that, in practice, the number of stars occulted every day 
by each asteroid is much smaller than one, so to decide whether 
an occultation event take place or not, we generated a random 
number uniformly distributed between zero and one, and we 
took as positive events those cases where the random number 
was found to be greater or equal to the number of occulted 
stars. If the test was positive, the epoch of the occultation event 
was generated randomly inside the time step interval. We re- 
peated this procedure every day over a period of 5 years, in 
order to consider different apparitions for each MBA, thus in- 
cluding different sky patches visited by each body. 

We assumed our observer to be situated at the geocenter. 
This choice is not a loss of generality, since the overall statistics 
will be the same for any given site on the surface of the Earth. 
However, we divided by a factor of 2 the number of occulta- 
tions obtained, to take into account the half celestial sphere vis- 
ible in the presence of a theoretical, geometric flat horizon. One 
should note that this correction could be different (and proba- 
bly larger) in practice, due to the difficulty of observing events 
at low sky elevations. Moreover, to calculate asteroid positions 
we solved the unperturbed two-body problem: planetary per- 
turbations represent a level of refinement not required for our 
purposes. As for asteroid diameters, the same values assigned 
as described in Sect. 12.21 were adopted. Finally we computed, 
for each of the simulated events, the maximum duration of the 
occultation r c using the object angular diameter and apparent 
motion at the occultation epoch, the asteroid's apparent magni- 
tude m a , and the magnitude drop during the occultation. 

The obtained results, of course, have a statistical value. 
They are presented in Fig. [2] where we show the average num- 
ber of occultation events that occur for a single site and a sin- 
gle object, during one year. It is clear that including faint stars 
strongly enhances the occultation rate. In fact, we find that any 
object larger than 20-30 km occults a 20' /! magnitude star at 
least once per year. 

However, the theoretical occultation rate that we have com- 
puted so far represents a theoretical value based on geometry 
only. It does not take into account the "observability" of each 
event, that is the capability of a given instrument to reveal the 
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Fig. 2. Number of occultations per object during one year, as 
a function of the asteroid diameter (thin lines). These statistics 
are valid for a single site and are computed for different values 
of the limiting magnitude as indicated by the legend. The thick 
curve (scale on the right-hand side y axis) represents the num- 
ber of asteroids for each bin in the simulated population. The 
product of this last number with each rate/object curve provides 
the total number of occultations per year for a given size range. 

temporary flux drop produced by the occultation. Furthermore, 
since the scientific goal that interests us is the impact of oc- 
cultation observations on asteroid size determination, we re- 
stricted our study to events from which the the size of an aster- 
oid can be determined within a given uncertainty level. 

We investigate the observability of the events by compar- 
ing portable equipment, similar to those currently used by most 
observers in present days, to a completely different situation, 
i.e. that of a fixed 1-m class telescope equipped with a cooled 
CCD camera. With the larger instrument, much fainter stars 
(and much smaller magnitude drops) can be measured, extend- 
ing the investigation of occultation events toward smaller aster- 
oids. As we will see, this will in turn allow the exploitation of 
improved asteroid ephemeris and stellar astrometry. 

4. The measurement of an occultation: limits 
imposed by the instrument 

The observation of an occultation event consists of the mea- 
surement of the temporal flux variation, by means of rapid pho- 
tometry, of a single source composed by the unresolved images 
of the star and the asteroid. A flux drop occurs during the oc- 
cultation when only the asteroid hiding the star is observable. 
The duration of the occultation is typically of a few seconds for 
large (50-100 km) asteroids, and the most favorable events, of 
course, involve a bright star and a faint asteroid. In this case, the 
flux drop can reach several magnitudes. We will use the sym- 
bols m a , m s , and m sa to indicate respectively the magnitude of 
the asteroid, the magnitude of the star and the magnitude of the 
unresolved source star + asteroid immediately before and af- 
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ter the occultation. For the corresponding measured fluxes we 
adopt the notation F(m a ), F(m s ), and F(m sa ), respectively. 

To obtain useful size estimates of asteroids from the ob- 
servation of an occultation event, the use of performing equip- 
ment, capable of exposure times considerably shorter than the 
maximum event duration, and sensitive enough to allow accu- 
rate photometric measurements of the flux drops is required. 
Still, a single observation of an occultation will not yield, in 
general, a direct measurement of the object's size. This is be- 
cause the measured occultation chord does not necessarily cor- 
respond to the largest projected size of the asteroid occulting 
the star. Moreover, in general, the three-dimensional shape of 
the occulting asteroid and its orientation should be known for a 
meaningful interpretation of the results obtained from occulta- 
tion measurements. 

As for the definition of the signal to noise ration (S/N) 
of the photometric measurements of the occultation event, 
we adopt the ratio between the amplitude of the flux drop 
(F(m as ) - F(m a )) and the photometric uncertainty associated 
with the flux during (cr> (a )) and outside (cr> (flJ )) the occultation: 

F(m„„) — F(m„) 
S/N = ) — - (2) 

y F(a) T u F(as) 

For a given instrument setup and for each event, the avail- 
able flux and the imaging equipment characteristics determine 
the corresponding photometric uncertainty. 

After having chosen an instrument set (see below) we thus 
computed - for each event - the minimum exposure time t m 
providing S/N~3. Assuming that the absolute timing is avail- 
able for the CCD frames with negligible uncertainty, we take 
t m as the estimate of the occultation timing uncertainty, i.e. 
cr T = t m . Since the maximum transit duration is proportional 
to the object diameter D, the ratio cr = ct t /t (where t„ rep- 
resent the maximum occultation duration, on the center line) 
will represent the relative accuracy expected for the diameter 
determination. 

One should note that, according to the S/N definition above, 
an exposure time yielding S/N~3 will always be found, but in 
some cases it will be comparable to - or significantly longer 
than - the occultation duration. However, we will consider as 
"non-interesting" those events having a relative uncertainty on 
the final diameter determination cr D > 40%. By adopting this 
constraint, "unobservable" events (i.e. those requiring t m < r ) 
are discarded. We then classified the remaining events as func- 
tion of (To- 

Of course, discarded events could still be of scientific value, 
for example for the search of asteroid satellites, however in this 
study we do not address in detail this otherwise interesting sub- 
ject. 

The predictions of occultation events that are currently 
distributed to observers include target stars not fainter than 
V=12-13, thus reducing the occultation rate at about 0.1 
events/objet/year even for the largest bodies. As a further se- 
lection, requirements on the minimum occultation duration and 
magnitude drop mean that to only the ~1000 asteroids larger 
than ~40 km are taken into account. We ran a first simula- 
tion considering the performance of a 20 cm, f/10 telescope 



equipped with a rapid, uncooled video camera (30 frames 
per second, exposure 1/30 sec) with rather high readout noise 
(15 e /pixel). We considered a CCD operating at 290 °K, at a 
pixel scale ~0.9 arcsec/pixel, with a seeing FWHM around 1 
arcsec. Sky background was assumed to be V=12 mag/arcsec 2 . 
With the high CCD frame rate imposed, the dominating noise 
source is not dark current, but the read-out noise. The sensor 
temperature is thus of secondary importance. By using aper- 
ture photometry with an optimized aperture size, we obtain the 
results shown in Fig. [3] The integral of the occultation rate for 
objects d > 40 km yields no more than a few events observable 
at high accuracy each year for a single site. This appears to be 
fully consistent with the present observation statistics. 
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Fig. 3. Number of observable occultations per object during 
one year, as a function of the asteroid diameter, for a 20 cm 
telescope equipped with a cheap fast-readout camera. These 
statistics are valid for a single site and are computed for dif- 
ferent values of the maximum relative uncertainty on size, cr D , 
as indicated by the legend. 

The simulation has also been run considering a 1-m tele- 
scope coupled to a low-noise, photon-counting CCD camera. 
These cameras, based on an electron multiplier in the read- 
out register, operate at very low temperatures (~ 180 °K) and 
presently represent the state-of-the-art technology for astro- 
nomical use. Considering a larger pixel (22 pmi) and a faster 
focal ratio (f/6), the image scale (0.76 arcsec/pixel) is not far 
from that of the 20 cm telescope equipped as above. Read-out 
noise is considered to be (1 e~/pixel); all other conditions are 
not varied relatively to the previous numerical experiment. The 
increase in camera performance and photon flux provides an 
improvement in the number of events at a given size range, by 
a factor of 1 to 2 orders of magnitude (Fig. |4j». 

In the following section, we use the statistics obtained with 
these two extreme instrument setups as a reference, and we in- 
clude in our analysis the contribution of uncertainties on both 
asteroids and star positions. 
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Number of occultations/object observable at different accuracies 
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Fig. 4. Number of observable occultations per object during 
one year, as a function of the asteroid diameter, for aim tele- 
scope equipped with a cooled, low-noise, fast-readout camera. 
The histogram is computed, as in Fig. [3] for different values of 
a- D . 

5. Role of astrometric accuracy 

Up to now we have discussed the statistics of occultation events 
and their observability as if we were able to know asteroid and 
star positions with infinite accuracy. In reality, they are affected 
by astrometric uncertainties. As a consequence, in principle, an 
observer has to monitor a larger number of candidate events in 
order to reach the count of the positive occultations predicted 
by the theoretical rate. The number of events that need to be 
monitored can be estimated by assuming as "impact section" 
for each asteroid the apparent size given by its diameter plus 
the width of the ephemeris uncertainty at the occultation epoch. 
Adding an uncertainty on the star position will further increase 
the "impact section". 

To take this effect into account, for each event we computed 
an estimate of the total uncertainty on the occultation path as: 

(Tpch = Jo%xj + o~l, where cr CEU is assumed to be equal to 
the CEU for the asteroid discussed above, while <x t represents 
the uncertainty on the star position. The prediction lists usu- 
ally diffused to the amateur community include the stars of the 
catalogues FK6 (I+III), Hipparcos, Thyco 2, and UCAC2 for 
stars of magnitude 10<V<12.5 (Dunham & Herald 2005I) . For 
stars with V<10 the UCAC2 positions, if available, are used 
to replace those given in the Thyco 2. A reasonable estimate 
of star position uncertainties is of about 1 mas for the FK6, 
10 mas for the Hipparcos, 70 mas for Tycho 2, and 20 mas for 
the UCAC2, for sources in the range of magnitudes V~10-14, 
whereas the uncertainty is of about 40 mas outside this range 
(Zacharias et al. 2004). For each simulated event we thus at- 
tribute to the occulted star an uncertainty typical of Thyco-2 or 
UCAC2 stars, with equal probability (this roughly reflects the 
current distribution of catalogue exploitation for predictions, 
see lDunham & Herald 20051 . 



In a second step, we compute the ratio p = 20/cr pat i t , (we 
recall here that 26 is the apparent asteroid diameter) averaged 
this over all the events occurring in a given diameter range. 
This can be considered to be a "prediction efficiency factor". 
If we assume that the uncertainties on asteroid ephemeris and 
star positions are Gaussian distributed, then cr /M ,/, represent the 
standard deviation of the total uncertainty on the occultation 
path, and p = 1 implies a 1-sigma confidence level (68%) of 
actually being able to observe the occultation event. This could 
certainly be considered as an event deserving an observation, 
with high success probability. In reality, the CEU can follow 
a more complex distribution, and the a- pat h has to be taken as 
a rough estimate of dispersion rather than a true Gaussian de- 
viation. Thus p cannot be considered more than a reasonable 
guess of the success probability for each event. On the other 
hand, one could also note that l/p, for p < 1, is an estimate 
of the number of events that should be observed to witness one 
positive event, on average. For those events with p~l the oc- 
cultation prediction is accurate and, in those cases where we 
have p » 1 the occultation path is known so accurately that 
it could be even possible to determine the position of the chord 
relative to the object barycenter. 

In the case of the after-Gaia statistics, different assump- 
tions are needed for both CEU and <x„. Since the dependence 
of the ephemeris accuracy on the asteroid size should pre- 
serve the general features shown in Fig. Q] we simply as- 
sume an improvement factor of 100 in the ephemeris uncer- 
tainty of each object, i.e. a post-Gaia value of the CEU re- 
duced by a factor of 100 with respect to the present one. This 
is a reasonable assumption in agreement with different analyt- 
ical estimates and simulations of the improvements in aster- 
oid orbital elements that Gaia will allow ( Virt anen et al. 2 004 
IMuinonen et al. 20041 IHestroffer & Berthier 20051 We also 
assume that the end-of-mission accuracy on star positions (in 
mas) can be represented by: 

crjmas) = 4 x 10~ 2 e °.44(y-20) + i.35 > y > 12 ( 3 ) 

where V is the star magnitude. A value cr* = 17 pas is obtained 
atV=150. In this case cr, saturates to the value reached at V=12 
for brighter stars. 

The result of the computation of p is given in Fig.[5]for the 
present and post-Gaia asteroid ephemerides and star position 
uncertainties. It is clear that the present occultation predictabil- 
ity appears reasonable (p>0.1) for the largest Main Belt aster- 
oids only (D>100 km), but in the "post-Gaia" era an equivalent 
accuracy can be reached down to ~15 km. 

6. Discussion and perspectives 

Let us consider a 10-20 km asteroid as a representative exam- 
ple. By multiplying the number of objects in the model at the 
corresponding bin (10 4 ) with the current number of occultation 
per object per year for a small instrument (Fig. [3), one finds 
that about 0.1 occultations/year occur at a 5-10% level of size 

4 This accuracy is obtained from simulations of the on-board de- 
tection ( Arenou & Lim 2003 ) and from the expected mission perfor- 
mances. 
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Fig. 5. Prediction efficiency factor for "present day" and "post 
Gaia" accuracies on asteroid ephemeris and star positions. 

accuracy for a single site. Furthermore, prediction efficiency 
(Fig- El) is so low (~0.002) that several hundred events should 
be observed (on average) to obtain one success. It is thus sim- 
ply impossible, with this equipment, to obtain accurate asteroid 
occultation diameters in this size range. 

On the other hand, for asteroids with D ^100 km the pre- 
dictability is rather good (~0.2) and bright star occultations 
(V<12), observable with the same equipment, occur ~0.01 
times per object per year (at high accuracy). Monitoring sev- 
eral objects by several sites will be sufficient to obtain, over 
a period of years, a sample of observed positive events. This 
is compatible with a lower limit at ~100 km for the currently 
available occultation-derived diameters. 

One should note that observing occultations with larger di- 
ameters can improve the accuracy of photometry, and increase 
the number of candidate events by including smaller flux drop. 
However, with present days asteroid ephemerides uncertainty, 
the predictability of the occultation events remains poor, thus 
preventing efficient observations at small sizes. 

Much more interesting is the case of the post-Gaia sce- 
nario, in which asteroid orbits are improved by a factor ~100 
and stellar astrometry makes a revolutionary step forward for 
star brightness down to V ~ 20. In fact, in this case the pre- 
diction efficiency remains p= 1 (or larger) down to a diameter 
of 20 km. This implies not only a perfect predictability of the 
events themselves, but also an uncertainty lower than the width 
of the occultation shadow on Earth. In other words, observers 
of occultations of asteroid with D >20 km will have some in- 
dications of the position of the chord relative to the projected 
shape of the object. This represents a major change in the way 
asteroid occultations can be approached, and could make them 
an even more powerful instrument for asteroid studies. 

In Fig. [6] we provide the distribution of the occultation du- 
rations and star magnitudes for events involving asteroids in 
the 15-20 km range. Comfortably long occultations and rela- 
tively bright stars are still present even toward the small end of 
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Fig. 6. Distributions of the occulted star V magnitude (upper 
panel) and of the duration of the events (lower panel) for all oc- 
cultations involving asteroids in the diameter range 15-25 km. 

the asteroid size distribution for predictable occultations. Also, 
we computed that ~70% of the events have a flux drop larger 
than 1 magnitude. However, in order to exploit the full potential 
of asteroid occultations in the post-GAIA era, more powerful 
instruments must be devoted to this activity. Figure @] shows 
that - by multiplying the predicted rate by the corresponding 
number of objects - 20 to 50 occultations/year are theoreti- 
cally observable with the highest success probabilities for each 
10-20 km object, by using aim telescope, with diameter un- 
certainties <10%. All known objects in this size range will lie 
within the CEU/0 ~ 1 limit if a two-orders of magnitude or- 
bit improvement is reached after Gaia (Fig. [TJ. This result has 
impressive implications, since it means that, in principle, all 
objects larger than that limit could be measured by monitor- 
ing occultations. Furthermore, the relatively large number of 
events could allow a single site to collect good statistics for 
each object, thus allowing an accurate diameter determination. 
In practice, the large number of asteroids involved (10 4 ) im- 
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plies that a few hundred events should be monitored each day. 
Of course, this is a huge effort for a single instrument or team, 
but it remains feasible if spread over a few years or if several 
sites with similar instruments are cooperating. 

The number of events at smaller asteroid diameters will still 
be interesting, although completeness will be out of reach. A 
selection of objects with the lowest ephemerides uncertainty 
will represent a good sample for extending size distributions 
toward km-sized asteroids. 

Gaia will not only provide accurate astrometry, but also 
photometry and spectroscopy. Experiments of lightcurve in- 
version from spare data (ICellino e t al. 2006 ) show that shapes 
and rotation vectors could be retrieved for potentially all the 
main belt asteroids from Gaia observations. Earth-based sur- 
veys such as Pan-STARRS could provide a significant con- 
tribution, although photometric accuracy will be much lower. 
However, the availability of shapes and physical ephemerides 
will allow extraction of a maximum of information from occul- 
tation studies, since the orientation of the asteroid as projected 
on the sky at the occultation epoch will be known. 

Several asteroid satellites will also be discovered and their 
relative obits determined, thus enriching even more the poten- 
tial of further occultation events. 

In conclusion, all the statistical arguments presented in this 
study indicate that, when Gaia observations will be available, 
a small set of dedicated 1 m telescopes will be able to provide 
- within a few years - a complete census of asteroid diame- 
ters down to ~10 km in size. The equipment foreseen to be 
used for future asteroid occultation studies such as fast, low 
noise CCD cameras, represent the state-of-the-art present-day 
technology, but it is reasonable to predict that in ~10 years 
from now (the after-Gaia era) even more optimized and per- 
forming devices will be available. Today, low noise CCD cam- 
eras are beginning to be used in professional occultation equip- 
ment dSouza et al. 2 006 ) and we can expect that they will en- 
counter a wide diffusion in the future. Considering the constant 
progress made in parallel by instruments accessible to the non- 
professional community, it is also highly likely that the role of 
dedicated amateur astronomers will continue to be a precious 
source of scientific information for asteroid studies. 
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